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ABSTRACT 

This thesis investigates the performance degradation resulting from 
multitone interference on orthogonal, noncoherent, frequency-hopped, M-ary 
frequency-shift keyed receivers (FH/MFSK) where the effect of thermal and other 
wideband noise is not neglected. The multiple, equal power jamming tones are 
assumed to correspond to some or all of the possible FH M-ary orthogonal 
signaling tones. Furthermore, the channel is modeled as a Ricean fading channel, 
a possibility precluded when thermal noise is neglected. Both the signaling tones 
and the multiple interference tones are assumed to be affected by channel fading. 

Both band and independent band multitone interference are considered. 

Performance is evaluated by obtaining a union boimd on the probability of bit 
error, and receiver performance is compared with exact results for band multitone 
interference of a noncoherent FH/MFSK receiver under comparable 
circumstances. Except for the case of Rayleigh fading of the signal, the union 
bound is very tight for those cases that can be compared with exact results. The 
advantages of the union bound approach are twofold. First, the union bound 
approach yields a solution that is far less computationally intensive than that 
obtained with the exact approach. Second, the union bound approach allows 
numerical results to be obtained for interference conditions that are not amenable 
to exact analysis, such as independent multitone interference of FH/MFSK. 
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I. INTRODUCTION 



A. SCOPE OF THESIS 

The goal of this thesis is to determine the electronic counter-counter measures 
(ECCM) potential of a non-coherent frequency-hopped M-ary frequency-shift keying 
(FH/MFSK) communication system under conditions of worst case hostile electronic 
counter measures (ECM) and fading channels. The performance of the FH/MFSK 
receiver will be investigated for frequency-nonselective, slowly fading Ricean channels 
and hostile multitone interference. 

The receiver considered uses a quadratic detector in each M-ary branch of the 
receiver for signal detection. Both band and independent multitone interference are 
considered and both jammer and signal are assumed to be affected by the fading channel, 
where it is also assumed that the channel fading need not necessarily affect the signaling 
tones and the interference tones in the same way. The cases of Rayleigh, Ricean, and no 
fading are all considered. 

Performance is evaluated by obtaining a union bound on the probability of bit 
error. The results are compared with the exact results, examined in [1], to determine the 
accuracy of the imion bound approach. Then the union bound is used to obtain probability 
of bit error for situations where the exact approach is either very difficult or even 
impossible to obtain. 



B. IMPORTANCE OF A SPREAD SPECTRUM COMMUNICATION SYSTEM 

The benefit from investigating the performance of a spread spectrum 
communication system lies in the extent to which spread spectrum communications are 
currently used. Their application to any military communication system aids in achieving 
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secure, low probability of intercept, and/or low probability of detection information 
exchange between friendly units as well as reducing the ability of hostile elements to jam 
friendly communications. Other applications of spread spectrum are: 

1) The use of either frequency-hopping or direct sequence spread spectrum in 
multiple access communications (as in cellular phones where vehicles communicate with 
a central station) in which a number of independent users are required to share a common 
channel without an external synchronization mechanism. 

2) The application of spread spectrum techniques in wireless local area networks 
(LANs). Wireless LANs are very important for the establishment of a network for mobile 
forces that circulates among the cooperating units various kinds of information, such as 
text messages, video, or audio. Wireless LANs are also very popular in commercial 
communications and are used for wireless communication of laptop computers with a 
control station. 

In any case, the potential of frequency-hopping spread spectrum to reject 
interference of either unintended or intended transmissions with the communicating 
signal is of crucial importance. The problem considered in this thesis is the evaluation of 
the performance of a frequency-hopped system in order to determine the system's 
capabilities and limitations. 

In chapters 2, 3, and 4, the union bound approach to performance analysis, a brief 
description of frequency-hopping, and the notation used in the thesis are discussed. In 
chapters 5 and 6, receiver performance with independent and band multitone jamming, 
respectively, are analyzed. In chapter 7 numerical results for both independent and band 
multitone jamming are presented. 
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II. IMPORTANCE OF UNION BOUND APPROACH 



The union bound approach simplifies the mathematical development by deriving 
an expression for the upper bound on the probability of symbol error that is always much 
simpler than the exact form. There are usually two approaches used to evaluate the 
probability of symbol error of a communications system. 

In the first approach, known as the indirect approach, the probability of symbol 
error is obtained using the relation 



where is the probability of not making a symbol error. Noting by X,,X2,...X^^the 
random variables representing the output of each branch of the receiver for a signal with 
M symbols, we have: 



The calculation of this probability is difficult, especially in the case where we deal not 
only with gaussian noise but also with jammer tones affecting the signal. This will 
become particularly obvious later on when we analyze receiver performance with 
independent multitone interference. In that case, the joint probability distribution 
functions required to solve the problem using the exact solution approach make it 
impossible to obtained a closed form expression for the probability of symbol error for a 
system using more than two symbols. 

However, if we use the imion bound, or direct approach, the solution works as 
follows, providing an upper boimd to the probability of error calculation: 



P = 1 -P 

S S 



( 1 ) 



P^ = 1 _ P^ = 1 _ Pr(X, > Xj n X, > X30...nX, > X^ll) . 



( 2 ) 
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P3 =Pr(X, <X2UX, <X3U...uX, <Xm| 1 ) < Pr(X,, <X2|l) + Pr(X, <X3|l) 
+...+ Pr(X, <X„| 1 ) 



( 3 ) 



We see that in place of joint probability distribution functions (pdfs) we now have the 
sum of individual pdfs. In fact, the individual pdfs can be grouped in groups of identical 
pdfs occurring with a certain multiplicity, so that the solution is a linear combination of a 
few pdfs. 

Summarizing, the union boimd approach provides a simple, fast method to 
investigate communication and jammer signal scenarios that are very difficult, if not 
impossible, to explore with the exact result approach as described in [ 1 ]. By using the 
union bound, we can increase our knowledge of commimication systems performance for 
a larger variety of jamming conditions. However, since this evaluation is a boimd, its 
accuracy, or tightness, needs to be tested so that we can comment on its effectiveness and 
limitations as a tool. In fact, in the following pages it will be shown that in some cases the 
union boimd approach is so accurate that it provides virtually identical results with the 
exact solution, while in some other cases the results obtained display an anomalous 
behavior. 
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III. BACKGROUND 



The main idea of frequency-hopping spread spectrum communications is to 
provide protection against a hostile jammer by increasing the bandwidth that the 
information signal occupies much more than the minimum required. By increasing the 
frequency range over which a system effectively operates, we force the jammer to spread 
its power over a wider frequency band and thus be less effective from a "per-frequency" 
signal corruption capability point of view. 

The main idea of frequency-hopping spread spectrum is based on the 
multiplication of the signal prior to transmission by an intermediate frequency, generated 
by a frequency synthesizer, of the form 

fj = f , -I- (i - 1) • Affl, where i=l ,2,...N (4) 

where N is the maximum number of possible frequency hop bins, Af^, is the separation 
between the carrier frequencies of adjacent bins, and 'i' changes pseudo-randomly every 
Tc seconds. By so doing, the entire spectrum of the signal transmitted, is shifted from its 
carrier frequency f,. to the new carrier frequencies: 

fc, =fc+fi (5) 

The result is that a new data modulated carrier that is shifted from one frequency to the 
next is created. From (5) we see that there are N different frequency hop bins, each of 
bandwidth Af^, . The value of i (and thus fc; ) is changed periodically according to some 
predetermined, but apparently random to a third party observer, noiselike spreading code 
called either the pseudo-random or the pseudo-noise (PN) sequence. 
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Generally, frequency-hopping is used with noncoherent modulation formats such 
as noncoherent MFSK. Such a signal is described by: 

s(t) = V2A, cos[27r(f5 + (m-l)A,)t + e,] (6) 

where m=l,2,...M, 0; is the i* symbol phase, and Af is the frequency separation between 
each of the M signaling tones. Multiplying this signal by a signal at the intermediate 
frequency f; = f, + (i - 1) • Afj^ , we get: 

s'(t) = 2 V 2 A, COS {H^s + (m - 1) Af ]t + 0j } cos(27ifjt) = 

= >/2A^ cos|2:r[fj -i-fj +(m-l)Af]t + 0i|-i- V2A^ cos|27t[fj -f^ -(m-l)Af.]t + 0j| 

( 7 ) 

The carrier frequencies of the first part are 

fi + f^ + (m - l)Af = f, + (i - 1)A4 + f^ + (m - l)Af (8) 

and is smallest for i=l and m=l; i.e., f; + f, . The carrier frequencies of the second part 
are 

fi - fs - (m - l)Af = fi + (i - l)Affl, - f3 - (m - l)Af (9) 

and is largest for i=N and m=l; i.e., 

f,+(N-l)Af„-f.. (10) 
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If the smallest frequency of the first part is greater than the largest frequency of the 
second part, then 

2f, >2f, +2B + (N-l)Affl, =^f, >f, +B + (N-1)A4 /2 (11) 

where B is the required guardband above and below the high and the low frequency 
signaling tones respectively. In this case, the signal s'(t) can be high-pass filtered to 
remove the frequency difference contribution, and the frequency-hopped signal becomes: 

V2A, COS {2n[f,-H+(m-l)A,]t + e,) (12) 

One obvious advantage is that a hostile jammer needs to know the pseudo-noise 
sequence of the communicator's transmitter in order to jam the specific hop bin where the 
transmitter operates at each time instant. To overcome a follower jammer, the signal must 
hop to a new carrier frequency before the jammer is able to complete its tracking process. 

Another advantage is due to the fact that the power spectral density of the 
frequency-hopped signal is identical to the that of the conventional signal in a specific 
hop bin. However, since the signal hops from bin to bin, and assuming that the 
probability that any bin is occupied is equal to (1/N), the average power spectral density 
is: 



PSD=^-2;s„,(f|f.=f.J. (13) 

i=l 

Hence, the signal power spectral density is lowered by a factor of N on average. 

Before proceeding we need to distinguish between the two main tone jammer 
methods used against communication systems. The first is band multitone interference 
and the second is independent multitone interference. In each case, we assume that the 
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multiple jamming tones are all of equal power. In band multitone interference, the 
multiple equal power jamming tones are distributed randomly across the entire frequency 
hop bandwidth, but the number of jamming tones placed within each jammed frequency 
hopped band is specific. In independent multitone interference, the multiple equal power 
jamming tones are distributed in a completely random fashion, meaning that the number 
of jammer tones within a frequency hop bin can vary from zero to a maximum equal to 
the total number of signal tones corresponding to each symbol the communicator can use. 
We also assume that each jamming tone corresponds exactly to one of the possible 
M X N signaling tones. From the communicator's viewpoint, this is the worst case. 

Before proceeding to the analysis of the receiver performance with multitone 
jamming, the notation used in the remainder of the thesis is discussed in the next chapter. 
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IV. FH/MFSK SYSTEM, CHANNEL MODEL, AND SYMBOLOGY USED 



The FH/MFSK communication system examined is assumed to have N non- 
overlapping FH bins. The bandwidth of each one is B Hz, where B is the bandwidth of 
the non-hopped MFSK signal. So the total frequency-hopped bandwidth is augmented to 
be greater than or equal to B-N. We also express by M the modulation order where 



M = 2' 



(14) 



and k is the number of bits per M-ary symbol. Since the number of signal tones within a 
hop bin is M and the total number of bins is N, the result is a total of M-N possible signal 
tone locations. 

The symbol duration T^ is related to the hop duration by: 



Te=KT^ 



(15) 



where K=l,2,3... . The symbol and hop rates are the reciprocal of the corresponding 
durations, so 



K= — 

' T3 



(16) 



and R = -^ , 



(17) 



respectively. Symbol rate is related to hop rate by 



9 



Rs=KR, (18) 

Using the frequency hopped signal (12), we note that the M symbols are assumed 
to be equally likely, both fi and fs are taken to be integer multiples of the symbol rate Rs, 
and the frequency difference between signal tones Af is taken to be an integer multiple of 
the symbol rate. Hence, 

A f = p • R^ (p integer). ( 1 9) 



To incorporate the idea of fading channels, we model the amplitude Ac as a Ricean 
random variable. This includes Rayleigh fading as a specific case. 

The transmitter and receiver are assumed to function in perfect synchronization, 
and the receiver after dehopping the signal uses a bank of M quadrature detectors where 
the integrator time constants are normalized to the symbol duration for notational 
convenience (Figure 1). 

The existing noise interference is expressed as A WON with a flat average power 
No 

spectral density defined as . The total power of the jammer is expressed as Pj^ , and 



the jammer delivers q equal power interfering tones such that the power of a single 
interfering tone is equal to 




( 20 ) 



The subscripts used in the description of the Ricean random variables are 'c' for 
the signal parameters and 'j' for the jammer. Thus, the pdf of Ac is 









•exp 



2a^ . 



•lo 



/ 2 2 ' 
a c -a c 






•u(aj 



( 21 ) 
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where u(«) is the unit step function, Iq(») is the modified Bessel function of the first kind 
and order zero, a is the power of the direct signal component, and 2a^c is the power of 
the diffuse signal component of the respective tone. The total average power of either a 
Ricean faded jammer or signal tone is equal to the mean squared value of the amplitude a^ 
or aj, respectively, and is 

P. =a^+2a^ (i = c,j). (22) 



The average power is assumed to remain constant from hop to hop. It is convenient to 
define the signal-to-noise power ratios of the direct signal components as 



Pi = 




(i = cj) 



(23) 



and the diffuse signal components as 



^i = 




(i = cj ) 



(24) 



where ci^„ is the noise power. 

For orthogonal signaling, the average probability of bit error is related to the 
average probability of symbol error by [2] 



_ M/2 

” "m-T 



(25) 
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and the average energy per bit is related to the average energy per symbol by 
E, =log2(M)-Eb =k-Eb. (26) 



12 



V. FH/MFSK RECEIVER PERFORMANCE WITH INDEPENDENT 

MULTITONE INTERFERENCE 



With independent multitone interference, the multiple, equal power jamming 
tones are distributed randomly across the entire frequency-hopped bandwidth. Two 
possible cases exist concerning the signal and jammer tones within a specific frequency- 
hop bin. In the first case there is no interference at all in the bin (the jammer tones are 
elsewhere), and the probability of error is the same as for a conventional MFSK signal in 
the presence of AWGN only. In the second case there is jammer interference in the bin. 
Since a bin is composed of M different signal locations and a potential jammer may 
place one or more of its tones to the aforementioned signal locations, this case is frirther 
divided into two subcases. In the first subcase the signal and jammer tones are on 
different branches, while in the second subcase one jammer tone coincides with a signal 
tone. This gives a total of three individual cases that are treated below. 

Case 1: 

In the following, we express by Ri, R 2 ,...Rm the random variables representing the outputs 
of branches 1,2,. ..M, respectively, of an MFSK receiver. Without loss of generality, we 
assume that the symbol T is transmitted. In order not to have an error, the output of 
branch one must be greater than the outputs of all other branches. The opposite situation 
expresses the error condition. As in equation (2), the probability of symbol error provided 
by the union bound is 

Pj = Pr(R, < Rj u R, < R 3 U...UR, < Rm|i) 

<Pr(R, <R 2 |l) + Pr(R, < R 3 |l)+...+ Pr(R, <R^|1) 

<(M-l) Pr(R, <R 2 | 1 ) ^ ^ 

and the last equation is true since the individual probabilities are equal. 
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Case 2a: 



Suppose that we have L jammer tones in the hop bin where the signal is. By Sj we express 
the existence of a jammer tone, and the subscript after j corresponds to the specific branch 
where the jammer tone is present. For example Sj^ means that a jammer tone is present in 
branch one. To simplify the discussion for this case, suppose that we have two jammer 
tones (L=2). Now 




= Pr(error with no interference in signal branch) 



(28) 



and Pr(R, < ) = Pr(R, < R 3 |l,Sj- ). 



(29) 



then (28) reduces to 



Pr(error|l, J < 2 Pr(R, < R^ |1, ) + (M - 3) • Pr(R, < R^ |1) 



(30) 



Generalizing for L jamming tones per hop bin, 1^<M-1, we get 



L 



Pr(error|l,£Sj,J<L-Pr(R, <R 3 |l,S^.J + (M-L-l)-Pr(R, <R^|1) 



(31) 




( 32 ) 
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As before, this can be simplified to 



Pr(error|l,Sj,,S 3 ^)< Pr(R, <RJl,Sj_,SjJ + (M-2) Pr(R, <R, 11 , 83 ,). 



(33) 



Generalizing, for L jamming tones per hop bin, 1<L<M-1 , we get 



Pr(error|l,Si, .Js,, ,, ) = (L - ])■ Pr(R, < Rjll.S;, ,S,^ ) + (M - L) Pr(R, < R„|l,Sj, ) 
1=2 

where it is understood that the summation is zero for L = 1 . 



(34) 



Now the probability of symbol error with at least one jamming tone per hop bin (L>1) is 
given by: 



Pr(error)= Pr(error|l. Z ) • Pr(l, X Sj,., ) + Prlerrorll.Sj, . Z Sj,., ) ' Prd, S , Z Sj,., ) 

1=1 1=1 1=2 1=2 

(35) 



L 

z 

1=1 



L 

z 

1=2 



L 

z 

1=2 



Since we have L jammer tones and M signaling tones per jammed frequency-hop bin. 



Pras„.zs*.,)— 



(36) 



and Pr(l,yS; )= 1 . 



1=1 



(37) 



So the probability of symbol error for L jamming tones per hop bin, L>1, is given by 
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( 38 ) 



P.(L)<i.[(L-l).Pr(R, <R;|l,Si,,sJ + (M-L).Pr(R, <R„|1,S„)]- 



+ [l-^j-[L Pr(R, <R,|l,sJ + (M-L-l)-Pr(R, <R„|l)] 



We now look for the probability of having L jamming tones per hop bin. This 
leads to the expression for the total probability of bit error. For independent multitone 
interference, and since we assume q jamming tones randomly interfere with a total of 
N-M possible signal tones, l<q<N-M, the probability that a frequency-hop bin contains 
one interfering tone is 



q 

N-M' 



(39) 



The probability that a frequency hop bin contains a second interfering tone is 



q-1 

N-M-1 ■ 



(40) 



The probability that a frequency hop does not contain an interfering tone is 



1 - 



q 

N-M 



(41) 



The probability that a frequency hop bin does not contain a second interference tone is 



1 - 



q-1 

N-M-1 ■ 



(42) 
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From (39) - (42), we derive the probability that a frequency hop bin contains no 
interference tones as 




q 

N-M-1 




— 




-^1 

N-M-i/ 



(43) 



Similarly, the probability that there is one jamming tone in a frequency hop bin is 



N-M-(M-l) 



(44) 

In the same maimer, we define the probability that there are two jamming tones in a 
frequency hop bin as 



M-N 

q 



1 - 



q-l 



N-M-1 

M-i/ 



1- 



q-1 



N-M-2 



■n 1- 



M-NifV N-M-i 



q q-i 

m-n’n-m-i 



1 - 



q-2 

N-M-1 






1 - 



q-2 
N-M-2^ 



, 1 - 



q-2 

N-M-(M-1)> 



=n 



i=0 



q-1 

M-N-i 




q-2 ^ 
M-N-J ■ 



(45) 



Combining (43), (44), and (45), we get the probability that L jamming tones are in a hop 
bin, where 1<L<M-1, as 
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Pr(L jamming tones per hop bin) = 



q-(L-l)-l 

MN-(L-l)-l_‘ 



q q-1 q-(L-l) 

m-n'm-n-i‘ 'm-n-(l-i) 

q-(L-l)-l ~ 




q-i 

M-N-i 




q-(L-l)-l^ 

M-N-i J 



( 46 ) 



The situation where L=M must be treated separately. Looking at the two 
individual products multiplied together in (46), we see that the first product has L terms, 
corresponding to the probability that L out of the M signals of the bin are jammed, while 

j 

the second product has the remaining M-L signals of the bin which are not jammed. Since 
when L=M all the signals are jammed, for that case we need only the first product term. 
So the probability of having L=M jamming tones in a hop bin is 



L-l 

Pr(L = M) = n 



i=0 



q-i 

M-N-i 



(47) 



Now, the total probability of symbol error for 0<L^ is given by 

Pr(error) = Pr(L = O) - ?^{\) + Pr(L = l) - ?^{l\L = l)+...+ Pr(L = M) - P 3 (l|L = m) 

(48) 



The maximum number of jamming tones we can have in a bin is M if M < q 

and q if q < M. This can be represented as min(q,M). By combining equations (46), 

(47), and (48), the probability of symbol error is given by the relation 
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fM-1 



n 1- 



\ 



. i=0 

I min(M-l,q) 



N-M-y 

L-l 



• (M-l)Pr(R, <R 2 ll) 



+ 



L=1 



n 



q-l 



UM-N-i 



i=L ^ 



M-N-i J 



Ps(L)-^ 






[ 

n 



q-l 



UM-N-i 



7 -Ps(M) 



( 49 ) 



The probability Ps(M) is obtained by letting L=M in (38). This yields 

Pj(M) = (M - 1) • Pr(Rj < R 2 |l.Sj_ ,Sj^ ) . Based on the above probabilistic description, it 

is apparent that we need to obtain the following probabilities for the MFSK receiver: 

Pr(R, <R,|1), Pr(R, <R 2 ll,S^, Pr(R, < R,|1,S^^). 

and Pr(R, < RjIhSj.^ ,Sj^ ). 

A. Pr(R, <R 2 |l) 

For noncoherent, orthogonal BFSK, Pr(R, < Rjll) is equivalent to the probability of 
bit error: 



2 



Pb(ac) = T-e 



(50) 



where 




(51) 



For noncoherent, orthogonal MFSK, the only thing that changes is that T^ is replaced by 
T : 

^ s* 
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^2 _N^ 

^ n “ rp 



( 52 ) 



To account for the existence of a Ricean fading channel we multiply equation (50) with 
(21), the Ricean pdf, and then integrate over all possible values of a^ . This yields 



00 00 2 

Pb = JPb(ac)-fA,(ac)da, = J-exp 



2a. 



•exp 



-a. 



2a 



c 



K- 



exp 











a% 


f 

1 


2 


CJ^n 



•exp 



a^c 

2a^c 



''a. -a ^ 
' c 



y 



•u(aj 



a. 






(53) 



Let X = a^c ^ dx = 2a^da^, and (53) becomes: 



Pb = 



1 


[a/1 


00 




f > 

1 1 




T 


f l~\ 
a,Vx 




2ct^c 


• |a^ exp 
0 


2 


an a ^ I 




’^0 


1 J 



da. 



(54) 



We then use the identity [3]: 



00 m-hn 

Jx 2 •e"'^ • J„(2pVx)dx 
0 



mip^ exp 




m+n+1 

a 




where 3"m(*) is the Laguerre polynomial and is defined [3]: 



(55) 
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fH 



m 



V a y 



= ^ , I J 

p = o p' '^^-p^ 



\ aj 



(56) 



In our case m=n=0 and the whole polynomial evaluates to unity. In place of "a" we have 



1 1 



^2 ^ _2 
Q n c / 



ja 

while in place of "P" we have — f . So, 

a. 



Pk = 



1 



4g c 
1 



•exp 






g^c 

2a^c 

g^c 

2a^c 



Jexp(- gx) • Jo(2p>/x)dx = 



o 2 2 \ 

2a cCT n 



2 2 
.a c +C7 



•exp 



a 



2 ^/^ O/T- 2^2 > 

2a a n 



.4a. 



2 2 
+a „y 



2 + 2 



-^•exp 
a c 

a^n 



g^c 



2c' 



1 

^2 2 

g c +a n 



2 + 2 



-^•exp 
C c 

g^n 



g^ 



2c c +g n 



.(57) 



Using (23) and (24), we have the definitions; 



g 



Pg = direct signal-to-noise-ratio = 

g n 



2c' 



= diffuse signal-to-noise-ratio = — — 

g n 



Applying these definitions to (57), we get 



Pr(R, <R 2 ll) = ^^-exp 



2 + ^c 



(58) 
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B.Pr(R, <R,|l,Sj,) 



This probability of error is the same as that derived in section 4. A except that a 
jammer tone is on the signal branch. Since signal and jammer tone exist on same branch, 
we can consider that they constitute a new kind of (altered) signal existing on the branch. 



In other words, the power of the signal and jammer tones add together but not 
necessarily in phase. We then integrate with respect to 0j after multiplying by its pdf to 

eliminate the dependence on 0. The phase angle 0j is assumed to be uniformly distributed 

in [0,2ti] . 

Without jamming, the mean values of the random variables representing the outputs of 
the in-phase and quadrature components of the first branch are given by; 



So the probabilistic analysis yields the same result with the only difference being that the 
"old" signal power a^c is replaced by 



a^c +a^j +2a^aj cos0j. 



(59) 



Xj = y/la^ cos(0) 



(60) 



X,,, =V2a,sin(0) 



(61) 



and in this case the output of the first branch is given by: 




( 62 ) 



This results in the probability of error: 
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(63) 



Pr(R, <R 2 |l) = ^exp 



■ 2N J 



With januning we have for the in-phase, quadrature components, and their squared sum 
as the final branch output the following relations, respectively. 



Xj_ = >/2a^ cos(0) + -s/2aj cos^0j) 


(64) 


sin(0) -i- y/la- sin^0j) 


(65) 


and X^i, +X^q, = 2 a^c +a^j +2a^aj cos(0-0j)j 


(66) 



We define 0 j = 0 - 0 j as the phase of the jammer with respect to the received signal to 
get the final form of equation (59). Considering the equations above, we can understand 
why we substitute the previous signal power a^c in Pr(R, < RjIfiSj,) by the composite 

power a^c + a^j -i- 2a^aj cos^0j j . Substituting (59) into (63), multiplying by 
fg. (0 j ) = j and integrating over 0 j , we get 



1 

Pr(R, <R 2 ll,Sj,) = — • Jexp 



a^c + a^ j + 2a^a j cos(0 j j 



2N. 



d0; 



(67) 



We then use the identity 
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( 68 ) 




z • cosi 



(4,-e)]d(t, 



to obtain 








= — • exp 



V(a'c+a^) 



2a „ 





( ■ai'l 


•lo 


c j 




\ o n y 



(69) 



However, this probability of error is conditional on both aj and a^ . To eliminate the 

dependence on a^, we multiply by its pdf (Ricean) and then integrate for all possible 
values of a^.. After multiplying (69) by (21) and integrating, we get: 



Pr(R, <R,|l,Sj,) = 



•exp 






exp 



-J 0 



2a^ 



(a'c) 



•exp 



a^c +g^c 
2a^c 






a -a 

c c 

. a^ 



da. 



1 






a^c 


20^“’’ 




• exp 


2a^c 



00 




0 



2 2 
O' c O’ n 

2a^ca^n 






(70) 



We now have the identity[3]: 
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( 71 ) 



« 1 
jx • exp[- ax^ ] • I„ (px) • J„(yx)dx = ^ • exp 



2 .,2 



P -Y 

4a 






The correspondence of the terms in (70) with the ones in the identity are: 



^ - 2 2 jp ~_2 ’ y — 



J«c 

ZG qG n ^ n ^ c 



So the probability of error, after eliminating the dependence on a,,, is 



Pr(R.<R,|l,sJ = j^^f;^.exp 



1 a^j +a 



2 ^ 



I 2 G c+o\J 



•lo 



a. 



,2 2 “ i 

c +a „ ‘ 



(72) 



To eliminate the dependence on ^ , we multiply by its pdf (Ricean) and then integrate for 
all possible values of a,. The Ricean pdf for aj is given by: 










(73) 



Multiplying equation (72) by (73), we get: 



2{g^c +g^) 




•exp 






1 


a^j 




a=i+a=A'| 


r: 


1 


G^c +g\J 


2 ■ 
j 


1 2a^j ]) 




( a. •a '^ 








•lo 




J J 

^2 










V 


G j j 
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2(g^c +cj^ 



n) 



a. 



c + CJ n 



2 -^jl-Io 



•exp 

f 



^ 1 CT^ja^j +c^ja^c +cJ^ca^j +a^jG^„ +a^jG^c +a^jCT^„^ 






(g^ +g\)g^j 



«j-^j 



1 


^j 


1^ 


1 


G^j +(g^ +G^„)^ 


2 G^j(G^+G\) 


• 2 exp 

^ j 




2 ' 


(g^+G^)g^j ^ 



•exp 



1 G^jtt^c +(of.^jG^c +a^jG^n) 



V 



(g^c +CJ^ 



n)cj'j 



J 



a. 



01 <t2 , -.2 

^CJ c n 



•ajJ-Io 



a • a-^ 

j j 

V cj j y 



(74) 



To integrate (74) with respect to 3j we again use identity (71).The correspondence of the 
terms in (74) with the equivalent ones of the identity are: 



1 G^j+(G^+G^n) „ CX. . ttj 

a = ~- . , 2 \ 2 > P= 2 2 > Y = J-^ 

2 (g^c +CJ^n)cJ^j <J c+<y n O 

For the argument of the Bessel function, we have: 



a. 



2a 



. ^2 

a c + CJ n 









G^j +(g^c +<y^ 



n) 






2 . ^2 \^2 
4- a n JCJ j 



2 2 2 

G j +G c +CT n 



So the final Bessel function has the form: 



ac«j 



^2 . „2 . „2 , 

va j +a c +CJ n/ 



(75) 
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Now, 



2a 



1 

j "i" 1 2 2 2 

1 Va j + CT c + C7 „y 

2 (a^+a^„)a^j 



(g^ +g\)g^j 

g^j +g^ +g\ 



So the terms appearing in front of the integral of (74) reduce to: 



,„2 ,„2 

G c^"CJ n 



■ c g n ) ' g j 



g^n 



2 g^j(cj^+g^„) 



g^n 



2(g^j +g^ +g^) 



(76) 



Finally, the argument of the exp(*) function that the integral of (74) yields is: 



^ «c V (: 



4a 



Vg^c+cj^n^ V g^y 



2 +(g^ +g^) 

_2 (g^+g^)g^j 



and combining this with the leading exponential functions we get the final exponential 
function to be of the form: 



\(J C +C7 n/ V OjJ 

g^j +(g^ +g^n) 

_2 {g^+g^n)g^ 



1 g^j-g^ +(g^jg^ +a^jg^) 

2 (g^+g\)g^j 



2 



(g\ +g^j) 

g^j +(y\ +g\ 



(77) 



Combining (75), (76), and (77) as required by (71) we get: 
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^2 

Pr(Ri<R,|l,Sj,)= . ; v -exp 

ZyG j + a c + cj nj 



2 2 
a c +a j 



2 a^j+a^c+o^nJ 



acCtj 



0 ^2 , ^2 , ^2 
VCJ j + CT c + CT nJ 



( 78 ) 



Substituting 




in (78), we get 
Pr(R,<R,|l,Sj,) = 



1 


Pc+Pj ^ 


•lo 




^ 1 


5 i+ 4 =+ 2 '“'’ 


1 ^j+^c+2j 


Ui+^c+ 2 j 



( 79 ) 



C. Pr(R, <R,|1,S,^) 

We begin by designating V, and Vj to be the random variables representing the 
outputs of branch 1 (where the signal resides) and branch 2 (where the jammer tone 
resides), respectively. Their pdfs are given by 




1 


f ^ 2 ^ 

Vj +2a^c 






9 ■ 

2g n 




o 

-f 






1 


f ^ 2^ 

V2+2a j 


T 


^c-yFl 


9 * exp 

2a^n 


2a^n 

\ ^ ^ 


' ^0 





( 80 ) 



( 81 ) 
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Multiplying (80) by (21) and integrating, we get 




a^ 



2c „ -c 



“• exp |-^_2 

2 , o _2 



a, -exp 



-a^ 



2c c 

c^n +2c c 



2c^ 
f 



•exp[ - -^_2 

n 



V 



2c^c^ 



\ <7 n y 




We then use the identity given in (71) where 



(82) 



a = 



C^n +2-C^c 

2c^„c^c 




Y = j 




Now we have 



V^tXc 

P-Y _ c^c^ _ 

2a C^n+2c^c C^n+2c^c 



(83) 



(84) 



2vi g^c 

P^~Y^ c\ c\ 2viC'’c +g'‘na^ 

4a 2(c^n+2c^c) 2(c^n +2c c)cT^nCr^c 



(85) 



Hence, using (83), (84), (85) in (71), we can evaluate (82) to obtain 



fv.(v,|l) = 



2(c^ +2c^) 



•exp - 



1 V, +2a^c 

2 c^n+2c^ 



-ac " 
_C^ +2c^ ; 



(86) 
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Similarly, 



fv.fv.lSiJ" 



2(CT^„+2CT^j) \ 2 CT^n+2o^jy 



1 Vj +2a^j ^ ^ f 



■h 



a^„ +2a^j 



m 



Using the random variable transformation R^i = V, , we get 






^ 1 f| "I" 2ct^c 



2 G \+2a cJ 



•I 



V2r, -a^ 

0l 2 , o_2 

\ Q n c y 



( 88 ) 



and similarly, 






1 rj +2a^j ^ 

2 a^+2a^jv 



V2t2 -a^ 



O’ n + 2c" 



V 



(89) 



Since R, and Rj are independent random variables, the joint pdf of R, and Rj is equal to 
the product of the individual pdfs: 



fR.Rj (^1 »’"2 j - f^^ (r, |l)- fR^ (tj |Sj^ j 



(90) 



( 00 00 00 

R, <R,|1,S, U fK_dr,dr, = Jf,,(r,|l). 

Or, 0 



ffR,(r 2 lSJ 2 )dr 2 dr, (91) 



For the inner integral. 
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( 92 ) 



00 00 

|fK.(r2lSi.)dr, =1^5775;^ exp - 

*■1 *“1 _ 



Fj + 2a^j 

2(a^ +2c^j) 



T 


r V2r2a, ^ 


* -^0 


^< T ^„ +2c^j J 



Let X = , , ^ = =>dx = -p===i=== 

-y/a „ +2 ct j -y/o „ +2(7 



X-yjc^n +2c^j 

CT^n +2a^j 



(93) 



Also to transform the lower limit of integration, when t 2 = r, , 



yfa\+2^j 



(94) 



Using (92), (93), and (94) to (91), we get 





f 1 




f 1 


00 

Pr(R, <R 2 |l.SjJ= J x-exp 

ri 


x^ c^„+2a^j 


■lo 


" 2 “ 2 


r -~—’ * X 

l^-y/a^n +2a^j j 


7cy^n+2cy^j 


^ J 







(95) 



Now compare the integral above with the Marcum’s Q-function [2], [4]: 



x^ +a^^ 



Q(a,P)= jx-exp[- — - — ^J-Io(ax)dx, 



and notice that if we let 



ct j r, 

a= , , ^ ■ , P = 



-y/a^ +2a^j ’ + 2<J^j 



(96) 
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Then (95) is given by 



f 

Q(a,P) = Q 

\ 



„ + 2a^j ’ -^a^n + 2a^j > 



( 97 ) 



Substituting (97) in (9 1), we get 



Pr(R,<RAsJ = ]-5;^.exp 



J c „ + 2a c 

f 

Q 






r I + 2a 
2(a^„ +2a^). 



V2r, 



r.a. 



Va^n + 2 a^y 



,-^a^n +2a^j ’-^a^ +2a' 

Consider the following identity [5]: 



(98) 



j / 



V ^2 ^2^ 



R. 



^•exp 
g I 



a^+R^ 

2a^i 



■I, 



^ a,R, ^ 

a^i y 



dR, = 



g^ 



g^i +a^2 



If we let 



1-Q 









2 ’t _2 



y^ya i+a 2 i+g 2 



a^ 



^2 , _2 

g 1 +g 2 



I '1'2 


I "2 


I a 2 I 


I a 1 


V 2 , 2 ’V 




l^V g I +g 2 V 


g 1 +g 2 j 



(99) 



g 2 =i/g^n +2a^j , a^=a-->l2, R, =r, , g, = -y/g^n +2a^ , ^i^=^f2■a^ 
then (99) can be used to evaluate (98) with the result 
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Pr(R,<R,|l,sJ = :^ 

+2p^j 

2(g^„ +g^ +g^j) 



g^n + 2g^ 



+ g^ +cf^j) 





r 1 


1 

1 

/O 


J 



a^c 



a 



.2 5 



an+CJc+C7j Van+CJc+CT^ 



j y 



r 


1 


1 O?o 


\ 


.\ 


|g^ +g^ +g^j ’) 


1 g^n +g^ 


+ g^j J 



( 100 ) 



Using (23) and (24), we get 



Pr(R, <R,11,SJ = 
. 1 + 



1 + ^c 
2 + ^c+^j 



1-Q 









V2 + ^,+4j’V2 + ^c+^jJ 



2 + ^c+^j \V2 + ^c+^j V2 + ^,+^J 



We can now use the identity [5]: 



Q(a,p) = l-Q(p,a) + exp 



a^+P^^ 



•Io(a,P) 



( 101 ) 



(102) 



in (101) to get 
Pr(R, <R 2 ll,Sj,) = Q 



V^Pj ■\/2Pc 



^2 + ^,+^. ’^2 + ^,+^. 
1 + ^, 



jy 



2 + 5.+5i 



•exp 



^ Pc+Pj ^ ^ 



V 



2 + ^.+4jy 



•I, 






2y^ 

v2 + ^c+U 



(103) 
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D.Pr(R,<R,|l,S,,,S,J 



This probability of error is obtained fromPr|^R, < R 2 |l,Sj J by altering the signal 
in branch one to include the jammer tone Sj^ . As before, a^c is replaced with 

a^c + a^j + 2a^aj cos^Gj) . We must then integrate twice to eliminate the dependence on 

the fading parameters. The phase angle 0 requires a third integration. However, the 
integrations cannot be evaluated analytically. The triple numerical integration is time 
consuming. As will be explained later, an effort is made to approximate 0 so as to avoid 
the third integration. 

The Pr|R, < R 2 |l,Sj^ j for a non-fading channel is given in [6] as 




< R2ll,Sj2) — 



1-Q 



_2 ’1-2 
Vi^n V^ny 




(104) 



Using (102), we get 



1-Q 



1 2 


1 2 ^ 


f 


2 1 7 ) 


/ 2 2 N 






a c 


a j 


= Q 




a j 


a c 


-exp 


a c + a j 


•lo 


a • a • 


■\l_2 M 


1 ^2 


J 


2 


1 2 


c 

rN 

b 

1 


2 




nj 




IV 


a n 1 


nJ 






^ a „ y 



(105) 



and substituting (105) into (104), we get 



p'(R'<r>|'-s-j=4J5’J5 



-2 exp 






V 



2a n 



^ n / 



( 106 ) 



Now, replacing a^c with a^c + a^j + 2a^aj cos(0) , we get 
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Pr(R, <R2ll,Sji,Sj2,a,,aj,e) = 



pr I 


[a^c +a^j +2a^aj cos(0) 


1 


^ a^c + 2a^ j + 2a,.aj cos(0)^ 




J 


^ 20" n ^ 






a^c + j + 2a^aj cos(9) 



V 






(107) 

The dependence on a^, aj , and 0 must be removed by multiplying (107) by the Ricean 

pdfs for the signal and jammer tones, respectively, and by the pdf for 0 and then 
integrating over a^, aj , and 0. As previously mentioned, this must be done numerically 

and can be a computationally intensive operation. In any event, the problem at hand has 
been "solved" in the sense expressed in [7]: 

In many cases we shall consider the solution found if we can express it in 
terms of one or more integrals, even if the integration can only be 
performed by expanding in series (or computing numerically). This is 
called 'reducing to quadratures', a phrase by which the mathematician 
symbolically washes his hands of the remainder of the task of finding the 
solution. 



This completes the analysis of the receiver performance for independent multitone 
jamming. The final form of the equations that must be evaluated are summarized in the 
appendix, equations B.6 through B.12. Many of these equations are also used for the 
analysis of band multitone jamming, which is addressed in the next chapter. 
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VI. FH/MFSK RECEIVER PERFORMANCE WITH BAND MULTITONE 

INTERFERENCE 

Band multitone jamming implies that the multiple, equal power jamming tones 
are distributed randomly across the entire frequency-hopped bandwidth with a specific 
number of jamming tones placed in a frequency-hop bin given that a particular hop bin is 
jammed. For worst case band multitone jamming, the worst probability of bit error is 
obtained when at most a single interfering tone exists at one of the M possible orthogonal 
signaling frequencies in a specific FH band. This section considers the general case of a 
jammer using L jamming tones per hop bin with 1^<M. We require the probability that 
the jammer places a set of L tones v/ithin the signal frequency-hop bin. This probability is 
derived by considering the following numerical example. Suppose that we have a total of 
10 jammer tones(q = 10), and two jamming tones per jammed bin (L = 2). Suppose that 
the number of hop bins is 30(N=30). The jammer tones are grouped into q/L = 5 sets, 
each set composed of two tones. Now, the probability that one of the five sets exists in 
one of the N = 30 hop bins is 5/30. In the general case, the probability that the hop bin is 



jammed is 




(108) 



The probability that the hop bin is not jammed is 



Pr(b in not jammed) = 1 



q 



(109) 



N-L 



From the total probability law, we find that the overall probability of error is: 
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Pr(error)= Pr(L, L > l) • P,(L, L > l) + Pr(L = O) • P,(L = O) 



( 110 ) 

The probability of symbol error when a hop bin is not jammed, Pj(L = O) , is the same as 
the probability of symbol error for noncoherent MFSK and no tone interference 

P,(L = 0)<(M-l)-Pr(R, <R2|1). (Ill) 

This non-jamming case must be considered in order to accoimt for every possible 
jamming scenario. 

The probability of symbol error when a hop bin is jammed is found in (38) to be 

(L>1): 



P,(L)£i [(L-l) Pr(R, <R;|l,Sj,,SjJ + (M-L).Pr(R, <R,|1,S,,)] + 

+ (l-^) [L Pr(R, <R,|l,sJ + (M-L-l)-Pr(R, <R,|l)] 

( 112 ) 

Hence, combining (110), (111) and (112), we obtain the probability of symbol error 
when (L>1) as 
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Ps(L) < — ^ 
^ ^ NL 



1 - 



'N-LJ 



(113) 



■g ■ [(L - 1) • Pr(R, < R, |1,S„ ,S J + (M - L) ■ Pr(R, < R, 11, Sj, )] + 
(l-^J-[L-Pr(R, <R,|l,S,J + (M-L-l).Pr(R, <R;|l)] 



(M-1) 



l_ 

2 + ^c 





The individual probabilities appearing in (113) are identical to the ones found in the last 
chapter, and are given in (58), (79), (103), and (107). The difference from the 
independent multitone interference is that L is fixed. 

This completes the analysis of the receiver performance for band multitone 
jamming. The final form of the equation that must be evaluated is given in the appendix 
as equation B. 13. In the next chapter, numerical results are obtained for receiver 
performance with both band and independent multitone jamming. 
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VII. NUMERICAL RESULTS 



The performance of the FH/MFSK receiver is investigated as a function of the 

P 

ratio of the information signal power to the total interference power, . Worst case 
performance is approximated by selecting q as the largest integer value between 1 and the 

Pj 

ratio . For non-worst-case performance, we select q=100. 

The non-fading, Ricean, or Rayleigh fading channel for the signal and/or the 
jammer tones is modeled by appropriately selecting the values of the corresponding 

a^c 

direct-to-diffuse power ratios, - - — and - ~ . A large ratio value implies little fading, 

2o^ c 2^ j 

a value of 10 implies Ricean fading, and a very small value implies Rayleigh fading. Due 
to programming necessities and in order to avoid divisions by zero, we use a ratio equal 
to lO”’ to simulate Rayleigh fading. For purposes of numerical computation, all curves 
are obtained with Eb/No = 13.35 dB and N=1000. The zero thermal noise curve is also 
plotted in each figure to validate the results since this provides an upper bound. For band 
multitone interference with a single interfering tone per hop and zero thermal noise, the 
probability of bit error is equal to [8] 



P =-3- 
N 
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fPj ^ 


TU 


— -Pc 
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1 q / 
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while for independent multitone interference with zero thermal noise, the probability of 
bit error for FH/BFSK is [9] 



2N ^ 
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A. BAND MULTITONE INTERFERENCE 



The performance of FH/BFSK with band multitone interference when the 
information signal is essentially unaffected by channel fading (a^c /2a^c = lOOO) with 

M=2 and M=8, respectively is shown in Figures 2 and 3. L successively takes the values 
1 and 2 in Figure 2 and the values 1, 4, and 8 in Figure 3. Performance is computed for 

P, 

both q = 100 and worst case q. When Pj > P^ , or equivalently ■— < -20dB and q=100, 

the worst probability of bit error is obtained for at most a single interfering tone at one of 
the M possible orthogonal signaling frequencies (L=l). For the region where P^ > Pj , the 

worst case performance is obtained when all signal tones are jammed (L=M). Hence, a 
jammer whose power Pj is guaranteed to exceed the signal power will do more harm with 
just one tone per bin. However, a jammer with Pj guaranteed to be less then the power of 
a signal tone, will do more harm by jamming all M tones of a bin. For worst case q, the 
worst performance is obtained when L=1 for all power ratios. 

The non-fading signal case for M=2 with all three possibilities of jammer tones 
fading is shown in Figure 4. The results are the same as those in [1]. When Pj > P^ or 



equivalently < -20dB and q=100, the difference in performance is trivial regardless 



Jt 



of whether or not the jammer tones are affected by the fading nature of the channel, even 
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when the jammer tones experience Rayleigh fading. On the other hand, when ^ Pj > or 

equivalently — > -20dB , the interference tones that experience Rayleigh fading yield 

the poorest performance. For the worst case q, the difference is trivial for the entire range 
of signal-to-jammer power ratios. Figure 4 is identical to the figure we get if we use the 
exact result equations. In fact, the individual probabilities of bit error used in the exact 
result, as they are described in [1], are identical to those obtained with the imion bound 
for the M=2 case. So M=2 will not be treated again by itself since it can convey no more 
information. It is important that for M=2 the simpler union bound equations provide the 
exact result. 

The non-fading signal case for M=8 is shown in Figures 5 and 6 for q=100 and 
worst case q respectively. The curves obtained for the exact result and the union bound 
are identical and practically indistinguishable from each other. The previous phenomenon 
with regard to the effect of fading interference tones deduced from Figure 4 is also 
obvious here, but it is even more pronounced, with the non-fading and Rayleigh fading 
curves differing by 5dB. 

The Ricean fading signal case for q=100 and worst case q is shown in Figures 7 
and 8 respectively, where the exact and approximate solution's curves obtained for M=8 
are compared. In Figure 7, we notice that the union bound curves are almost identical to 

P, 

the exact result ones and only in the region where > -lOdB do they appear to have a 

“h 

constant difference of 1.3 8-1 O''*. In Figure 8, the union bound curves behave as expected, 
being practically indistinguishable from the exact result ones for P, < Pj^ , while in the 

opposite case they appear to have a slight difference with a maximum value equal to 
1.38- 1 O'". The union bound curves again work as an upper bound on the probability of bit 
error. 

Summarizing the results up to this point, since the differences detected are trivial 
ones, for all practical purposes we infer that for M even larger than eight, we can easily 
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resort to the union bound equations to estimate the probability of bit error when the signal 
tones are either not affected at all by the fading channel or the signal tones experience 
Ricean fading. 

The exact result and the union bound curves for a Rayleigh fading channel 
affecting the signal tones and the three possible cases of fading conditions for the jammer 
tones are shown in Figure 9. The difference between the two groups of curves attains a 
maximum value equal to 0.0038. This difference between the exact and union bound 
curves is the greatest noticed up to this point. The union bound curves still act as a bound, 
being higher than the exact result ones, but the accuracy is questionable. This is examined 
further in the following figures where we determine if increasing the number of symbols 
M causes a decrease in the probability of error as should normally happen. 

The non-fading and Ricean fading signal curves in conjunction with every 
possible fading situation of the jammer tones is displayed in Figures 10 through 14. The 
value of M is sequentially changed, taking the values 2, 4, and 8. In these figures, no 
anomaly is detected, and the curves behave correctly in all cases with the probability of 
error decreasing as the number of symbols M increases. However, the next three figures. 
Figures 15, 16, and 17, display the Rayleigh fading signal case in conjunction with all 
three possible fading conditions of the jammer tones. In all cases the results obtained are 
incorrect. The union bound still acts as a bound, but the probability of error erroneously 
increases as M increases. 

Summarizing the union bound approach for band multitone interference, we 
see that for the non- fading and the Ricean fading signal case, the union bound approach 
gives very accurate results which are perfect in the first case and deteriorate some in the 
second though still providing acceptable accuracy. However, when the signal tones 
experience Rayleigh fading, the union bound approach should be avoided because the 
accuracy obtained is not sufficient. From a computational point of view, the union bound 
approach is far simpler than the exact approach. The exact approach equations include 
joint pdfs of two random variables expressed with double integrals. On the contrary, the 
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union bound equations are expressed as the sum of individual pdfs of one random 
variable and thus they are expressed in terms of one or more integrals. 



B. INDEPENDENT MULTITONE INTERFERENCE 

In independent multitone interference, the number of tones existing per frequency 
hop bin can vary from zero to M, where it equates to the number of signal tones. As 
mentioned before, the calculation of the probability of error here is time consuming due 
to the inability to further simplify the solution which contains three integrations. 

However, one integration is eliminated for the case where either the signal or the jammer 
tones are not affected by the fading channel, and two integrations can be eliminated for 
the case where both the signal and jammer tones are non-fading. 

The M=2 case for a non-fading signal and various jammer fading cases when we 
have independent multitone interference is shown in Figure 18. As with band multitone 
interference, the results obtained for M=2 are identical with the ones given by the exact 
solution. This is because the exact solution equations, described in [1], again coincide 
with the union bound for M=2, enabling us again to use the imion boimd with perfect 
accuracy. As in band multitone interference, we notice the performance degradation from 
the non-fading jammer case to the Rayleigh fading jammer case when > Pj . 

The non-fading signal case for M=8 is shown in Figure 19. It is obvious that the 
union bound approach is a valuable tool in this case since the exact result approach 
requires very complicated probabilistic equations with the complexity further increased as 
M increases. However, the union bound approach provides us with a concise expression 
for any value of M. The only means of comparing the curves obtained here with a pre- 
existing result is to get a feeling from the comparable ones obtained for band multitone 
interference case bearing in mind the result obtained in [1], that the probability of error 
for the band multitone interference is poorer than for independent multitone interference. 
In fact, comparing Figures 5 and 19, we see that they have the same shape, differing by a 
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factor of roughly eight. Also, comparing Figures 4 and 18, we see that they differ by a 
factor of roughly equal two. The difference between band and independent multitone 
interference is roughly equal to M. 

Figure 19 is characterized by the phenomenon previously described. For 
Pp > Pj a channel with no fading of the interference tones provides an advantage of up to 

5dB over one with Rayleigh fading. 

The Ricean fading signal case for M=2 with all three possibilities of fading of the 
jammer tones is displayed in Figure 20. As for non-fading signal tones, when q=100 we 
notice a difference in performance for the non-fading and Rayleigh fading jammer tones 
when P,. > Pj . That difference, however, is now less pronounced. For the worst case q, 

the difference between the respective curves is again trivial. Comparing this figure with 
Figure 3, we again notice that the probability of error is less by a factor of roughly two. 

The Ricean fading signal case for M=8 with all three possibilities of fading of the 
jammer tones is displayed in Figure 21. Comparing this figure with Figure 7, we notice 
that the probability of error difference is eight. This confirms the previous observation 
that the decrease in probability of error when we shift from band to independent 
multitone interference is roughly equal to M. 

The Rayleigh fading signal case for M=2 and M=8 'with non-fading interference 
for q = 100 and worst case q is shown in Figure 22. The obvious problem is the same as 
the one observed for band multitone interference. The accuracy achieved is not sufficient 
and predicts an increase in the probability of bit error for an increase in the number of 
symbols M. Thus, the union bound approach should be avoided when the signal 
experiences Rayleigh fading. 

The remaining figures. Figures 23, 24, 25, and 26, show the results of 
approximating 0 rather than integrating over all of its possible values. Remember that by 
0 we denote the difference between the phase of the signal and the phase of the jammer 
tones. This is modeled as a uniformly distributed random variable over [0, 2k]. As 
mentioned before, in the beginning of this section, for the case of more than one jammer 
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tone per frequency hop bin, we must perform three numerical integrations in order to 
calculate the probability of error: one for the amplitude of the signal tone fading, one for 
the amplitude of the jammer tone fading, and one for the angle 0. To this must be added 
the embedded calculation of the Marcum's Q-function that is also defined in integral 
form. It would be very convenient to find a way to eliminate one numerical integration. In 
fact, the integration over 0 was experimentally found to be the most time consuming one. 
After trying different angle values to get an accurate approximation, it is found that for 
7 ... 

0 = — • 71 the approximation is so good that the difference between the exact and 

approximated traces can barely be detected. The exact and approximated curves for a 
non-fading signal and various fading of the jammer tones for M=2 and M=8, respectively, 
are shown in Figures 23 and 24. The exact and approximated curves for a Ricean fading 
signal and various fading of the jammer tones for M=2 and M=8, respectively, are shown 
in Figures 25 and 26. This approximation results in a significant reduction in 
computational complexity and time. 
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VIII. CONCLUSION 



From a probabilistic performance point of view, the results obtained verify those 
of [1] that when Pj > the effect of fading jammer tones is very small when the signal 

tones are not affected by fading. When P^ > Pj performance is poorer when the 

interference tones experience fading. When the information signal experiences fading, the 
effect of fading interference tones is less. 

With regard to the imion bound approach, it is shown to be very accurate for non- 
fading or Ricean fading signal tones, while for Rayleigh fading signals its use should be 
avoided. 

The application of the union boimd for the Ricean fading signal in the case of 
independent multitone interference for M=2 and M=8, that was not treated in [1], shows 
that the probability of error advantage when the jammer uses independent multitone 
interference instead of band multitone interference is equal to a factor of M. Hence, when 
M = 8, independent multitone interference gives probabilities of error eight times less 
than does band multitone interference. 

In case of the independent multitone interference, the calculations involved can be 

7 

much simplified if instead of integrating over all possible values of 0 we let ® ^ ’ 

thus eliminating the most time-consuming integration. 

In summary, the union bound approach taken in this thesis not only simplifies the 
computational complexity involved in obtaining numerical results, but also allows us to 
analytically evaluate situations that are not amenable to solution using an exact approach. 
The principle drawback is the breakdown of the union bound approach when the signal 
experiences Rayleigh fading. 



49 





50 



2siii[(w, + {M - DAttijr] 

Figure 1 ; Noncoherent FH/MFSK Receiver 




L=1 

L=2 

L=1 worst case q 
“® L=2 worst case q 

Performance for zero thermal noise 



Figure 2: Performance of a FH/BFSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones for 
L=1 and L=2. 
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Pb 0 - 




“®" L=1 
L=4 
L=8 

“® L=1 worst case q 

L=4 worst case q 
L=8 worst case q 

Performance for zero thermal noise 



Figure 3: Performance of a FH/8FSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones for 
L=l, L=4, and L=8. 



52 




Pc/PjT(dB) 

"®” Union-bound result -jammer tone no fading 
Union-bound result - jammer tone Ricean fading 
■®" Union-bound result - jammer tone Rayleigh fading 
jammer tone no fading, worst case q 
jammer tone Ricean fading, worst case q 
jammer tone Rayleigh fading, worst case q 
“ Performance for zero thermal noise 



Figure 4: Performance of a FFl/BFSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones. The 
results are identical with those obtained using the exact result 
equations. 
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Exact result - jammer tone Ricean fading 
"®" Exact result - jammer tone Rayleigh fading 
Union-bound result - jammer tone no fading 
Union-bound result - j ammer tone Ricean fading 
Union-bound result - jammer tone Rayleigh fading 
Performance for zero thermal noise 

Figure 5: Performance of a FH/8FSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones, q = 
100, using both the exact result and the union bound equations. The 
curves obtained are identical and practically indistinguishable. 
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Pc/PjT(dB) 

Exact result - jammer tone no fading 
Exact result - jammer tone Ricean fading 
Exact result - jammer tone Rayleigh fading 
Union-bound result - jammer tone no fading 
Union-bound result - jammer tone Ricean fading 
Union-bound result - jammer tone Rayleigh fading 
Performance for zero thermal noise 



Figure 6; Performane of a FH/8FSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones, q„.<. 
using both the exact resxilt and the union bound equations. The 
curves obtained are identical and practically indistinguishable. 
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Exact result - jammer tone no fading 
Exact result - jammer tone Ricean fading 
Exact result - jammer tone Rayleigh fading 
Union-bound result - jammer tone no fading 
Union-bound result - jammer tone Ricean fading 
Union-bound result - jammer tone Rayleigh fading 

Figure 7: Performance of an 8FSK noncoherent receiver with Ricean fading 
of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 
q = 100, using both the exact result and the union bound equations. 
The union bound curves behave as expected, being a little higher 

Pc 

than the exact solution ones. For — > -lOdB we lose a little bit of 

pjj 

accuracy as the union-bound curves have the same shape with the 

exact result curves but they have a constant difference of 1.361 O' 

As in the M=2 case, there is no essential difference between no 
fading, Ricean fading and Rayleigh fading of the jamming tones, 
for Pj>Pc- 
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Pc/PjT(dB) 

Exact result - jammer tone no fading 
Exact result - jammer tone Ricean fading 
Exact result - jammer tone Rayleigh fading 
~o Union-boimd result - jammer tone no fading 
Union-bound result - jammer tone Ricean fading 
Union-bound result - jammer tone Rayleigh fading 
Performance for zero thermal noise 

Figure 8: Performance of an 8FSK noncoherent receiver with Ricean fading 
of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 
using both the exact result and the union bound equations. 

The union bound curves behave as expected, being practically 
indistinguishable from the exact result equations for Pc<Pj'p while 

for Pc>Pj j they appear a little bit higher, having a maximum 
difference of 1.3410 '’. 
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Pc/PjT(dB) 

Exact result - jammer tone no fading 
Exact result - jammer tone Ricean fading 
"®" Exact result - jammer tone Rayleigh fading 
“O' Union-bound result - jammer tone no fading 
“0- Union-bound result - jammer tone Ricean fading 
Union-bound result - jammer tone Rayleigh fading 

Figure 9: Performance of an 8FSK noncoherent receiver with Ricean fading 
of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 
q = 100, using both the exact result and the union bound equations. 
Again the fading behaviour of the jammer tone is irrelevant to the 
probability of error. The union bound curves have the same form as 
the ones from the exact result, but appear to display the biggest error 
than in any case treated up to is point, an error having a maximum 
equal to 0 . 038 . 
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M=4 

M=8 

"® M=2 worst case q 
M=4 worst case q 
M=8 worst case q 



Figure 10: Checking the Union Bound accuracy for essentially no fading of the 
information signal for band multitone interference with no fading of 
the jammer tones, for both q = 100 and , for M = 2, 4, and 8. No 
anomaly detected for the using the union bound and Pt decreases 
as M increases. 



59 




-e- M=2 

-0- m=4 

M=8 



“® M=2 worst case q 

M=4 worst case q 
-tJ- ivi=8 worst case q 

Figure 1 1 :Checking the Union Bound accuracy for essentially no fading of the 
information signal for band multitone interference with Rayleigh 
fading of the jammer tones, for both q = 1 00 and , for M = 2, 4, 
and 8. No anomaly detected for the using the union bound and 
Pb decreases as M increases. 
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M=2 

-0- M=4 

M=8 



”® M=2 worst case q 
M=4 worst case q 
M=8 worst case q 

Figure 12: Checking the Union Bound accuracy for Ricean fading of the 

information signal for band multitone interference with essentially 
no fading of the jammer tones, for both q = 100 and q,^,^ , for M = 2, 
4, and 8. No anomaly detected for the using the union bound and 
Pi, decreases as M increases. 



61 




M=2 

-0- m=4 

M=8 



M=2 worst case q 
M=4 worst case q 
M=8 worst case q 

Figure 13: Checking the Union Bound accuracy for Ricean fading of the 
information signal for band multitone interference with Ricean 
fading of the jammer tones, for both q = 100 and q,^^ , for M = 2, 
4, and 8. No anomaly detected for the using the union bound 
and Pb decreases as M increases. 
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M=4 
-B- M=8 

“® M=2 worst case q 

M=4 worst case q 
■B M=8 worst case q 

Figvire 14: Checking the Union Bound accuracy for Ricean fading of the 

information signal for band multitone interference with Rayleigh 
fading of the jammer tones, for both q = 100 and q,^,. , for M = 2, 
4, and 8. No anomaly detected for the Pausing the union bound 
and Pb decreases as M increases. 
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M=2 

M=8 

M=16 



M=2 worst case q 
"® M=8 worst case q 

“ M=1 6 worst case q 

Figtore 15 : Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
essentially no fading of the jammer tones, for both q = 100 and 
q^^j , for M = 2, 8, and 16. Anomaly detected!! ! For increasing M 
the probability of error increases. 
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M=8 

M=16 

“® M=2 worst case q 

M=8 worst case q 
M= 16 worst case q 

Figure 16 : Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
Ricean fading of the jammer tones, for both q = 100 and q^^, , 
for M = 2, 8, and 16. Anomaly detected!!! For increasing M the 
probability of error increases. 
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Pc/PjT(dB) 

M=2 
M=8 
-« M=16 

“® M=2 worst case q 

M=8 worst case q 
-D- m= 16 worst case q 

Figure 17: Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
Rayleigh fading of the jammer tones, for both q = 100 and q^^c , 
for M = 2, 8, and 16. Anomaly detected!!! For increasing M the 
probability of error increases. 
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o No fading jammer tone 
* Ricean F ading j ammer tone 
X Rayleigh Fading jammer tone 
+ No fading jammer tone worst case q 
-- Ricean Fadingjammer tone worst case q 
- - Rayleigh Fadingjammer tone worst case q 
. . . Performance for zero thermal noise 

Figure 1 8: Performance of a FH/BFSK noncoherent receiver with essentially 
no information signal fading for independent multitone interference 
with various conditions of fading of the multiple interference tones. 
The results are identical with the ones obtained using the exact 
result equations. 
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o No fading jammer tone 
* Ricean F ading j ammer tone 
X Rayleigh F ading j ammer tone 
+ No fading jammer tone worst case q 
“ Ricean Fading jammer tone worst case q 
- - Rayleigh Fading jammer tone worst case q 

Figure 19: Performance of a FH/8FSK non-coherent receiver with essentially 
no signal fading and various interference fading for independent 
multitone interference. 
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o No fading jammer tone 
* Ricean Fading jammer tone 
X Rayleigh F ading j ammer tone 
+ No fading jammer tone worst case q 
-- Ricean Fading jammer tone worst case q 
- - Rayleigh Fading jammer tone worst case q 
. . . Performance for zero thermal noise 



Figure 20: Performance of a FH/BFSK non-coherent receiver with Ricean 
signal fading and various interference fading for independent 
multitone interference. 
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o No fading jammer tone 
* Ricean F ading j ammer tone 
X Rayleigh F ading j ammer tone 
+ No fading jammer tone worst case q 

— Ricean Fading jammer tone worst case q 

- - Rayleigh Fading jammer tone worst case q 

Figure 21 : Performance of a FH/8FSK non-coherent receiver with Ricean 
signal fading and various interference fading for independent 
multitone interference. 
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Pc/PjT(dB) 



Figure 22: Performance of a FH/BFSK and FH/8FSK receiver for Rayleigh 
fading information signal and non-fading interference, for 
independent multitone interference. Anomaly detected, since 
for an increase in M we notice an increase to the probability of 
error. 
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Rayleigh Fading-exact 0 
Rayleigh Fading - 0 = 7 tu/ 1 0 
Non-fading, worst case q, exact 0 
Non-fading, worst case q, 0 = 7 ti/10 

Ravleieh fadine. worst case a. exact 0 
Rayleigh fading, worst case q, 0 = 7 ti/1 0 

Figure 23: A comparison of the union bound curves for a FH/BFSK 
noncoherent receiver's performance with essentially no 
information signal fading for independent multitone interference 
with either no fading or Rayleigh fading of the jammer tones, 
with q = 100 and worst case q, obtained by integrating 0 and also 

7 

by approximating 0 with — tx . The results obtained are identical, 

so the approximation can effectively be used simplifying the 
calculations. 



72 




Pc/PjT(dB) 

Non-fading - exact 0 
Non-fading, 6 = 7 tc/ 10 
Rayleigh fading - exact 0 
Rayleigh fading, 0 = 7 tc/ 10 
Non-fading, worst case q, exact 0 
Non-fading, worst case q, 0 = 7 tc/ 10 
Rayleigh fading, worst case q, exact 0 
Rayleigh fading, worst case q, 0 = 7 ti/10 

Figure 24: A comparison of the union bound curves for a FH/8FSK 
noncoherent receiver's performance with essentially no 
information signal fading for independent multitone interference 
with either no fading or Rayleigh fading of the jammer tones, 
with q = 100 and worst case q, obtained by integrating 0 and also 

7 

by approximating 0 with — ti . The results obtained are identical, 

so the approximation can effectively be used simplifying the 
calculations. 
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Non-fading, 0 = 7 ti/ 1O 
Raylei^ fading - exact 6 
Raylei^ fading, 0 = Tti^IO 
Non-fading worst case q, exact 0 
Non fading worst case q, 0 = 7 ti/ 1 0 
“*< Raylei^ fading worst case q, exact 0 
"® Raylei^ fading worst case q, 0 = 7 tc/10 



Figure 25: A comparison of the union bound curves for a FH/BFSK 

noncoherent receiver's performance with Ricean information 
signal fading for independent multitone interference with either 
no fading or Rayleigh fading of the jammer tones, with q = 100 
and worst case q, obtained by integrating 0 and also by 
7 

approximating 0 with— 7t . The results obtained are identical, so 

the approximation can effectively be used simplifying the 
calculations. 
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Pc/PjT(dB) 

Non-fading, exact 0 
Non-fading, 0 = 7 j:/ 10 
Rayleigh-fading, exact 0 
Rayleigh-fading, 0 = 7 ti/ 1 0 
Non-fading, worst case q, exact 0 
Non- fading, worst case q, 0 = 77t/10 
Rayleigh-fading, worst case q, 0 = 7jt/10 
■® Rayleigh-fading, worst case q, 0 = 7 ji/10 

Figure 26: A comparison of the Union Bound Curves for a FH/8FSK 

noncoherent receiver's performance with Ricean information 

signal fading for independent multitone interference with either 

no fading or Rayleigh fading of the jammer tones, with q = 100 

and worst case q, obtained by integrating 6 and also by 

7 ... 

approximating 0 with — 7t . The results obtained are identical, so 

the approximation can effectively be used simplifying the 
calculations. 
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APPENDIX . MATHEMATICAL IDENTITIES AND RELATIONS 



For convenience, all pertinent equations are repeated here. 

A. IDENTITIES 
1. Integrals 

• I „ (z) = exp[- • J „ ( jz) = (- j)" • J „ ( jz) 

• Io(z) = — • Jexp[±z-cos((l)-e)]d(l) 

Z7t Q 



m + n 



•lx 2 -J^(2pVx)dx = 



mip^ exp 
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a 
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,m + n + 1 
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(pZ 
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where is the Laguerre polynomial: 
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m (. 
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( 9 ^ 
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p = 0 
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U J 



CO j 

j X • exp[- ax ^ ] • I „ (px) • J „ (yx)dx = — • exp 
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.j f£il 

"V2aY 



(A.1) 



(A.2) 



(A.3) 



(A.4) 
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2. The Marcum's Q function and related integrals 



Q(a,p)= jxexp 



x^+a^ 



•Io(ax)d> 



(A.5) 



• Q(a,0) = l 



(A.6) 



>2\ 



Q(0,p) = expl -y 
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Q(a,p) = l-Q(p,a) + exp 
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•Io(a,P) 



(A.8) 
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(A.9) 



B. RELATIONS 



1. Union bound on the probability of bit error 

Pj = Pr(R, < Rj w R, < Rj u. . .uR, < Rm |l) 

<Pr(R, <R 2 |l) + Pr(R, < R 3 ll)+...+ Pr(R, <Rm| 1) 
<(M-l)-Pr(R, <R2|1) 

where by Rj, Rjv-Rm express the random variables representing the outputs of 
branches respectively, of an MFSK receiver. 



80 



2. pdf of a signal's or jammer's amplitude modeled as a Ricean random 
variable 





•exp 



2a'i J 










(B.2) 



where u(*) is the unit step function, Io(*) is the modified Bessel function of the first kind 
and order zero, a^c is the power of the direct signal component, and 2a^c is the power of 
the diffuse signal component of the respective tone. 



3. pdf of a signal's or jammer's amplitude modeled as a Rayleigh random 
variable 








(B.3) 



4. Signal-to-noise power ratio of the direct signal components 



O n 



(i = c,j ) 



(B.4) 



5. Signal-to-noise power ratio of the diffuse signal components 
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(B.5) 
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6. Pr(R, <R 2 |l) 
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9. Pr(R,<R,|l,Sj^,Sj^,a„a.,e) 
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(B.7) 



(B.8) 



(B.9) 



82 



10. Pr(R,<R,|l,Sj,,S,J 
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11. Probability of symbol error for L jamming tones per hop bin, 1<L<M 



P,(L)£i-[(L-l).Pr(R, <R,|l,Sj,,SjJ + (M-L).Pr(R, <R„|l,Si,)] + 
- ■ [l ■ Pr(R, < R, |1,S J + (M - L - 1) ■ Pr(R, < R„ |l)] 



(B.ll) 



12. Probability of symbol error for independent multitone interference for L 
jamming tones per hop bin, 0<L<M. 
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where the probability Pj(M)is obtained by letting L=M in Ps(L). This yields 
P,(M) = (M - 1) • Pr(R, < R^lhSj. ,S J . 
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So 
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13. Probability of symbol error for band multitone interference for L 
jamming tones per hop bin, L>1 
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14. Relation between probability of bit error and probability of symbol error 
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• Pj , where M is the number of symbols 
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